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Abstract 

On the basis of earlier investigations on homeoidally striated Mac 
Laurin spheroids and Jacobi ellipsoids (Caimmi and Marmo 2005; 
Caimmi 2006a, 2007), different sequences of configurations are de- 
fined and represented on the ellipticity-rotation plane, (Oex^). The 
rotation parameter, x^, is defined as the ratio, E rot /E res , of kinetic 
energy related to the mean tangential equatorial velocity component, 
M(v^) 2 /2, to kinetic energy related to tangential equatorial compo- 
nent velocity dispersion, and residual motions, M(a^ uw + 
o"3 3 )/2. Without loss of generality (above a threshold in ellipticity 
values), the analysis is restricted to systems with isotropic stress ten- 
sor, which may be considered as adjoint configurations to any assigned 
homeoidally striated density profile with anisotropic stress tensor, dif- 
ferent angular momentum, and equal remaining parameters. The de- 
scription of configurations on the (Oex^) plane is extended on two re- 
spects, namely (a) from equilibrium to nonequilibrium figures, where 
the virial equations hold with additional kinetic energy, and (b) from 
real to imaginary rotation, where the effect is elongating instead of 
flattening, with respect to the rotation axis. An application is made 
to a subsample (N = 16) of elliptical galaxies extracted from richer 
samples (N = 25, N = 48) of early type galaxies investigated within 
the SAURON project (Cappellari et al. 2006, 2007). Sample objects 
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are idealized as homeoidally striated MacLaurin spheroids and Jacobi 
ellipsoids, and their position on the (OeXv) plane is inferred from ob- 
servations following a procedure outlined in an earlier paper (Caimmi 
2009). The position of related adjoint configurations with isotropic 
stress tensor is also determined. With a single exception (NGC 3379), 
slow rotators are characterized by low ellipticities (0 < e < 0.2), 
low anisotropy parameters (0 < 5 < 0.15), and low rotation pa- 
rameters (0 < Xv < 0.15), while fast rotators show large ellipticities 
(0.2 < e < 0.65), large anisotropy parameters (0.15 < 5 < 0.35), and 
large rotation parameters (0.15 < x\ < 0.5). An alternative kinematic 
classification with respect to earlier attempts (Emsellem et al. 2007) 
needs richer samples for providing additional support to the above 
mentioned results. A possible interpretation of slow rotators as non- 
rotating at all and elongated due to negative anisotropy parameters, 
instead of flattened due to positive anisotropy parameters, is exploited. 
Finally, the elliptical side of the Hubble morphological sequence is in- 
terpreted as a sequence of equilibrium (adjoint) configurations where 
the ellipticity is an increasing function of the rotation parameter, slow 
rotators correspond to early classes (E0-E2 in the oblate limit and 
E— 2-EO in the prolate limit) and fast rotators to late classes (E3-E6). 
In this view, boundaries are rotationally distorted regardless of angu- 
lar momentum and stress tensor, where rotation has to be intended as 
due to additional kinetic energy of tangential equatorial velocity com- 
ponents, with respect to spherical configurations with isotropic stress 
tensor. 

keywords - galaxies: kinematics and dynamics - galaxies: elliptical. 

1 Introduction 

According to their original classification (Hubble 1926), elliptical galaxies are 
subdivided into eight classes, designated E0, El, E7, the numerical index 
being the integer nearest to lOe, where e is the ellipticity of related galaxies 
as projected on the sky, and < e < 0.7 is inferred from observation^. 
Hubble's original scheme has generally been considered satisfactory, but the 
physical explanation of the distorted boundary has been changed in time. 

The classical interpretation starts from the evidence that (i) all celestial 
bodies are known to be in rotation, and (ii) the symmetry of figure shown by 
elliptical galaxies is precisely such as rotation might be expected to produce; 

1 In modern classifications, galaxies within the class E7 are considered as lenticulars 
instead of ellipticals (for further details refer to e.g., Caimmi 2006b), but it is unrelevant 
to the aim of the current investigation. For this reason, Hubble's original classification 
shall be considered in the following. 
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which, in turn, suggests (iii) an inquiry as to how far the observed figures of 
elliptical galaxies can be explained as the figures assumed by masses rotat- 
ing under their own gravitation. In this view, boundaries are distorted by 
systematic rotation around the minor axis i.e. angular momentum parallel 
to the rotation axis. For further details refer to the parent textbook (Jeans 
1929, Chap. XIII, §§299-302) and to recent attempts (e.g., Caimmi 2006b). 

About fourty years ago, observations begun to yield increasing evidence 
that (giant) elliptical galaxies cannot be sustained by systematic rotation 
(e.g., Bertola and Capaccioli 1975; Binney 1976; Illingworth 1977, 1981; 
Schechter and Gunn 1979). Accordingly, (giant) elliptical galaxies were con- 
ceived as systems with triplanar symmetry and ellipsoidal boundaries set 
up by specific anisotropic stress tensor (Binney 1976, 1978, 1980) and a 
negligible contribution from angular momentum. Owing to high-resolution 
simulations, the same holds also for (nonbaryonic) dark matter haloes host- 
ing galaxies and clusters of galaxies (e.g., Hoeft et al. 2004; Rasia et al. 
2004; Bailin and Steinmetz 2004). In this view, boundaries are distorted by 
systematic rotation around the minor axis and/or anisotropic stress tensor. 

It is worth mentioning that an anisotropic stress tensor does not imply 
orbital anisotropy and vice versa. With regard to e.g., nonrotating systems, 
and leaving aside instabilities, orbits could be only radial or only circular, and 
the shape maintain spherical or flat. On the other hand, anisotropic stress 
tensors yield distorted boundaries (e.g., Binney 1976, 2005, the last hereafter 
quoted as B05). The extent to which angular momentum and anisotropic 
stress tensor are effective in determining the shape of a system, may be 
quantified by the ratio, x\i of kinetic energy due to figure rotation to kinetic 
energy due to remaining motions. It is, in turn, related to a squared velocity 
ratio, (x^)obs, inferred from observations. For assigned density profiles, the 
dependence of the rotation parameter on the shape may be determined using 
the tensor virial theorem. 

In the special case of homeoidally striated ellipsoids, where the isopycnic 
(i.e. constant density) surfaces are similar and similarly placed (Roberts 
1962), the rotation parameter, Xv, has been explicitly expressed as a function 
of the ellipticity, e, restricting to oblate shapes (Binney 1976). Configurations 
with assigned anisotropic stress tensor are represented on the (Oex v ) plane 
as a family of curves which branch off from the origin, (e,Xv) = (0,0), and 
increase for increasing ellipticity. The curve related to the isotropic stress 
tensor makes an upper limit, while their counterparts characterized by an 
anisotropic stress tensor (cr n = a 2 2 > 033) lie below (Binney 1976). 

The tidal action from an embedding, nonbaryonic dark halo makes the 
shape of the embedded elliptical galaxy, closer to its own. In other words, 
a less flattened dark halo makes a less flattened embedded elliptical galaxy, 
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and vice versa. Accordingly, a family of curves depending on two parameters, 
namely the fractional mass, m = M halo /M ega i, and the axis ratio of the dark 
halo (supposed to be axisymmetric), ehaio, can fit to the data on the (0e%„) 
plane even if the stress tensor is isotropic (Caimmi 1992). 

The rotation parameter, Xv, is independent of the density profile for 
homeoidally striated ellipsoids (Roberts 1962; Binney 1976) which, in gen- 
eral, make only a first approximation to equilibrium configurations (e.g., 
Vandervoort 1980b; Vandervoort and Welty 1981; Lai et al. 1993). With 
regard to equilibrium configurations e.g., polytropes (Vandervoort 1980a), 
the rotation parameter, Xv, depends on the density profile (Caimmi 1980). 
Accordingly, a family of curves depending on a single parameter, namely the 
polytropic index, n, can fit to the data on the (Oex v ) plane even if the stress 
tensor is isotropic (Vandervoort 1980a; Caimmi 1983). 

Sysytematic rotation and anisotropic stress tensor are effective to the 
same extent in flattening or elongating a spherical shape, provided imaginary 
rotation (around the major axis) is considered in the latter case (Caimmi 
1996b, 2007, 2008; the last quoted hereafter as C07, C08, respectively). The 
key concept is that the distribution function is independent on the sign of 
tangential velocity components, and the whole set of possible configurations 
is characterized by an equal amount of both kinetic and potential energy 
(Lynden-Bell 1960, 1962; Meza 2002). 

In other words, clockwise and counterclockwise tangential velocity com- 
ponents are indistinguishable to this respect, and no change in shape occurs 
for any configuration between the limiting cases: (i) equal amount of clock- 
wise and counterclockwise tangential equatorial velocity components, which 
implies a null mean value; (ii) only clockwise or counterclockwise tangential 
equatorial velocity components, which implies a maximum or minimum (ac- 
cording to the sign) mean value. In particular, a single configuration exists 
with isotropic stress tensor, to be conceived as "adjoint" to its generic coun- 
terpart with anisotropic stress tensor. The above considerations hold above 
a threshold in ellipticity (e > e a dj) which, on the other hand, lies below 
the threshold for elliptical galaxies (e > e e u > e ad j). Accordingly, elliptical 
galaxies always admit adjoints configurations with isotropic stress tensor. 

The tensor virial theorem provides a rigorous global link between angular 
momentum, anisotropic stress tensor and shape, but the quantities appearing 
therein cannot easily be deduced from observations. A classical approxima- 
tion to the rotation parameter, Xv, is (Xv)obs = V max /a , where V max is the 
peack rotation velocity and cr a centrally averaged velocity dispersion, both 
projected on the line of sight (e.g., Illingworth 1977, 1981; Schechter and 
Gunn 1979). With the advent of integral-field spectroscopy (Bacon et al. 
2002; de Zeeuw et al. 2002; Kelz et al. 2003) a rigorous connection can be 
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established between the tensor virial theorem and observations, provided the 
former is reformulated in terms of sky-averages. In particular, the rotation 
parameter inferred from the data is (x^)obs =< vf\ >^ 2 / < of\ >^ 2 , where v\\ 
is the mean velocity projected on the line of sight, a 2 is the related variance, 
and < v 2 >, < a 2 >, are sky-averages with respect to the mass (B05). 

In general, the expression of the rotation parameter depends on what 
is intended as rotation energy. Usually, the kinetic energy is decomposed 
into contributions from ordered and random motions (e.g., B05). In the 
current attempt, the kinetic energy shall be decomposed into contributions 
from cylindrical rotation, M(v^) 2 /2, and tangential equatorial component 
velocity dispersion, plus residual motions, M{p 2 ww + cr| 3 )/2; the 

rotation parameter, xl, shall be defined as the ratio of the first to the sum 
of the remaining two above mentioned terms. 

In this view, an extension to nonequilibrium configurations and imagi- 
nary rotation appears quite natural, in the light of a general theory where 
boudaries are distorted by angular momentum and/or anisotropic stress ten- 
sor (C07, C08) and triaxial configurations due to the occurrence of bifurca- 
tion points are also considered (Caimmi 1996a,b, 2006a,b; CM05). In dealing 
with imaginary tangential equatorial velocity components, the rotation pa- 
rameter, Xv, h as to be used instead of Xvi where positive and negative values 
are related to real and imaginary angular momentum, respectively. 

The present investigation is mainly devoted to the following points: (i) 
representation of nonequilibrium figures on the (Oexl) plane, restricted to 
homeoidally striated MacLaurin spheroids and Jacobi ellipsoids and edge- 
on orientations with the major axis perpendicular to the line of sight; (ii) 
location on the (Oexl) plane of elliptical galaxies with assigned inclination 
angle and anisotropy parameter, restricted to axisymmetric shapes and clas- 
sified as fast and slow rotators (Cappellari et al. 2007; hereafter quoted as 
SX); (iii) extent to which nonequilibrium figures fit to elliptical galaxies; 
(iv) interpretation of slow rotators as intrinsically prolate, nonrotating bod- 
ies, and related location on the (Oexl) plane; (v) representation of adjoint 
configurations in imaginary rotation with isotropic stress tensor. 

The work is organized as follows. The rotation parameter, xl> is expressed 
as a function of the intrinsic meridional ellipticity, e, for equilibrium and 
nonequilibrium figures with isotropic stress tensor, in Section [2j The position 
of elliptical galaxies on the (Oexl) plane, with regard to a restricted sample 
for which fiducial values of the inclination angle and the anisotropy parameter 
can be assigned, together with their adjoint configurations where the stress 
tensor is isotropic, is determined in Section [3j A comparison with earlier 
attempts is made in Section HI The possibility that slow rotators are in fact 
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prolate, nonrotating bodies, is exploited in Sectional where an interpretation 
of the elliptical side of the Hubble sequence is also proposed. The conclusion 
is drawn in Section El 



2 Nonequilibrium figures 

Given a mass distribution with assigned rotation axis, x%, the kinetic energy 
is usually decomposed into contributions from ordered and random motion 
(e.g., B05), where the former arises from mean velocity components within 
any infinitesimal volume element, and the latter from related variances which 
appear in the expression of the stress tensor. In this view, the contribution 
of streaming motion to velocity dispersion (along a selected direction) is not 
included in the stress tensor. 

A different approach shall be exploited in the current attempt. More 
specifically, the kinetic energy shall be conceived as the sum of two contribu- 
tions: one, related to either the mass-weighted tangential equatorial velocity 
component, M(v^) 2 /2, or the moment-of-inertia-weigthted angular velocity, 
Is(Q) 2 /2, and one other to the contribution of the related velocity dispersion 
plus the remaining velocity components, radial equatorial, (v%,), and radial 
polar, (v 2 ). The above contributions shall be hereafter quoted as rotation 
kinetic energy, E TOt , and residual kinetic energy, E ies , respectively. In this 
view, the contribution of streaming motion to velocity dispersion (along a 
selected direction) is included in the stress tensor. 

It is worth mentioning that the mass-weighted tangential equatorial ve- 
locity component and the moment-of-inertia-weigthted angular velocity yield 
different kinetic energy values (for a formal discussion refer to Appendix [Al. 
As the current investigation deals with linear velocities, only the former shall 
be considered here. 

If the tensor virial theorem holds for only time-averaged quantities, the 
related mass distribution can no longer be considered in dynamical or hydro- 
static equilibrium, but still in virial equilibrium (C07), which makes a special 
kind of nonequilibrium figures. In the following, "nonequilibrium" has to be 
intended as "far from dynamical or hydrostatic equilibrium but still in virial 
equilibrium". For nonequilibrium figures, the tensor virial equations must 
be related to an adjoint configuration with equal density profile and rotation 
energy, but different residual-energy tensor components, (E Tes ) pp = ( pp E Tes . 



The result is (C07): 



2(-Erot)gq -|- 2(^qgE TCS -\- (-Ep t)gq , 
2(33-Eres + (-Spot) 33 = ; 



9 = 1,2 
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where (E pot ) pq is the potential-energy tensor and the coefficients, ( pp , may be 
understood as generalized anisotropy parameters (CM05, Caimmi 2006a,b; 
C07). The combination of Eqs. ( 11 al) and ( llbl) yields: 



2(E 
2E 



.Jqq — 7^"(-^pot)33 + (E 



,33 



potjqcj 







1,2 ; 



C: 



33 



"(3 



pot J 33 i 



(2a) 
(2b) 



and the potential-energy tensor for homeoidally striated Jacobi ellipsoids 
reads (CM05): 



pot )pq 



GM 2 



^Bel\E se \)pq , 



p,q= 1,2,3 ; 



seljptj 



-Ssel — X/ e s2£s3^4s 



(3a) 
(3b) 



where 5 pq is the Kronecker symbol, G the constant of gravitation, z/ se i a 
profile factor (i.e. depending only on the density profile), a p are semiaxes, 
a p /a q axis ratios, A p shape factors (i.e. depending only on the axis 



-pq 



ratios). 

The substitution of Eqs. (j3J), ( 1421) . (j4"3l . into (TSJ) yields after some algebra: 



CM 
ai 



2,2,2 
<^ + + ^33 
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(,33 
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«1 C33 



C = Cn + C22 + C33 j 

C E TCS / E,. es 
Cpp {E res ) pp / E Tes Q 



pp 



( E Tes ) pp / _E r , 

(-^res)pp 



(4a) 

(4b) 
(4c) 

(4d) 
(4e) 



where the anisotropy parameters, ( pp , are related to equilibrium figures (CM05; 
C07; C08). 

The rotation parameter, xl — E Iot /E res , by use of Eqs. (HI) and ( 1411 takes 
the explicit expression: 
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which, for assigned values of the anisotropy parameters, depends only on the 
axis ratios. In an earlier attempt (CM05) a different definition of rotation 
parameter has been used, in terms of a mean squared velocity, < v^ ot >, 
instead of related squared mean velocity, < t> ro t > 2 , which implies differ- 
ent anisotropy parameters for equal shapes and rotation parameters, in the 
two formulations. The current definition is more strictly connected with ob- 
servations i.e. mean velocities and velocity dispersions. The factor within 
brackets in Eq. (l5al) corresponds to equilibrium figures = 1, ( pp = ( pp ) 
while the virial index, ( = E res /E res , is a measure of the departure from 
equilibrium (0 < ( < +00). 

For homeoidally striated Jacobi ellipsoids, the anisotropy parameters, £11, 
(22, are related to the diagonal components of the moment-of-inertia tensor, 
hi, I22, and a necessary and sufficient condition for the stress tensor to be 
isotropic, is £33 = 1/3. In addition, bifurcation points are independent of the 
angular momentum and the stress tensor, and may be related to the adjoint 
configuration where the squared angular momentum attains a maximum or 
minimum value (according to the sign) and the stress tensor is isotropic. 
This is why the effect of an anisotropic stress tensor is equivalent to an 
additional real or imaginary rotation, inducing flattening (on the equatorial 
plane) or elongation (on the rotation axis), respectively. For further details 
refer to the parent paper (C07) and Appendix [B] More specifically, real or 
imaginary rotation take place when the factor within brackets in Eq. floHI) is 
positive or negative, respectively. 

Given a nonequilibrium figure with assigned anisotropy parameters, ( pp , 
the particularization of Eq. ( l5al) to the adjoint configuration with isotropic 
stress tensor (£ 33 = 1/3) yields: 



(X?j)iso C 



1 -Bscl 



3(3 



selj33 



(6) 



which can be plotted as a function of the meridional ellipticity, e = 1 — 631, 
where e > implies real rotation, e = no rotation, and e < imaginary 
rotation. Cases within the range, 1/2 < ( < 2, are plotted in Fig.fU where 
the full curve corresponds to equilibrium figures (( = 1) and the upper on the 
first quadrant (( = 1.1, 1.2, 2.0) and the lower on the first quadrant (( = 
1/1.1, 1/1.2, 1/2.0) dotted curves to nonequilibrium figures. Points lying 
on the first and the third quadrant represent systems in real and imaginary 
rotation, respectively. The origin is the locus of nonrotating systems. The 
loci of different bifurcation points are represented by different vertical lines, 
namely: from axisymmetric to triaxial configurations (dotted) ; from triaxial 
(left) and axisymmetric (right) to pear-shaped configurations (dashed); and, 
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in addition: onset of dynamical instability in axisymmetric configurations 
(dot-dashed). The first and the third quadrant can completely be filled with 
curves expressed by Eq. (EJ), where the limiting case, £ = 0, corresponds to flat 
configurations with no centrifugal support, and the limiting case, ( — > +00, 
corresponds to nonflat configurations with no pressure support. 

In the general case of anisotropic stress tensor, conformly to Eq. (I5ap . the 
trend on the (Oexl) plane is similar but the curves are shifted on the right 
(nonrotating flattened configurations) or on the left (nonrotating elongated 
configurations), according if <j\ x = a\ 2 > O33 or o\ x = o\ 2 < <t| 3 , respectively. 
In any case, the adjoint equilibrium configurations with isotropic stress tensor 
lie on the related sequence (full curve in Fig. [I]) with unchanged axes, which 
implies a vertical shift of a selected point on the (Oexl) plane, until the 
adjoint configuration is attained. 

3 Comparison with observations 

Aiming in comparing model predictions to data from observations, the posi- 
tion of elliptical galaxies on the (Oexl) plane shall be determined. The sam- 
ple used (Caimmi and Valentinuzzi 2008; Caimmi 2009b; hereafter quoted as 
CV08 and C09, respectively) is extracted from richer samples of early-type 
galaxies investigated within the SAURON project (SIV, SX). The follow- 
ing parameters may directly be deduced or derived from observations: the 
effective (half-light) radius, R e ; the stellar mass within the effective radius, 
M e = M(R e ) = M/2; the luminosity-weighted average ellipticity, on a plane 
perpendicular to the line of sight, < e± >, within either an isophote en- 
closing an area, A = ttR^, or the largest isophote fully contained within 
the SAURON field, whichever is smaller; the luminosity-weighted squared 
mean velocity component, parallel to the line of sight, < t>| >, within ei- 
ther an ellipse of area, A, ellipticity, e±, and related position angle, or the 
largest similar ellipse fully contained within the SAURON field, whichever is 
smaller; the luminosity-weighted squared velocity dispersion, parallel to the 
line of sight, < erjj >, within either an ellipse of area, A, ellipticity, e±, and 
related position angle, or the largest similar ellipse fully contained within the 
SAURON field, whichever is smaller. For further details refer to the parent 
papers (SIV; SX) and an earlier attempt (B05). 

Two additional parameters can be inferred by fitting the data with dy- 
namic models. More specifically, the inclination angle, i, is deduced from 
the best fitting two-integral Jeans model (SIV), and the anisotropy param- 
eter, 5, is determined from the solution of the dynamic models, supposed 
to be axisymmetric (SX). It is worth remembering that fast rotators show 
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evidence of large anisotropy and axial symmetry, while slow rotators appear 
to be nearly isotropic and moderately triaxial (SX). 
The anisotropy parameter (e.g., B05): 

6=l-?f = l-?f ; (7a) 

*11 ^22 

o"n = o"22 > 0-33 ; (7b) 



may be related to the generalized anisotropy parameters, ( pp , via Eqs. (lid 
lej), in the special case of equilibrium figures, as: 

^ _ ^ C33 t C33 1 — 3^33 3£n — 1 3^22 — 1 /c 



C11 C22 1 — C33 C11 C22 

which is restricted to homeoidally striated MacLaurin spheroids. 

For representing sample objects on the (Oexl) plane, intrinsic values of 
velocities and ellipticities with regard to edge-on configurations have to be 
used, instead of projected values along the line of sight. The corrections for 
edge-on (edo) configurations are (Binney and Tremaine 1987, Chap. 4, §4.3): 

i-4 = 1 ~ < 2 ± >2 ; (9) 



1 - 631 = 1 " 



1- < e_L > (2- < e_L >) 
sm 1 



1/2 

; (10) 



for the meridional axis ratio and ellipticity, where < ej_ > and < e± > are the 
luminosity averaged axis ratio and ellipticity related to an inclination angle, 
i, between the symmetry axis and the line of sight (i = 90° for edge-on 
configurations), and: 

r ~2 1 [<^!>]obs 

< vi > cdo = L ; (11) 

11 sm 1 

1 2 1 [<<i>]obs 

< 0\\ > edo = Z j — ! 12 

11 1 — COS % 

under the assumption of axisymmetric [a\ = 02) configurations with axisym- 
metric (an = 022) stress tensor (SX). 

The intrinsic squared mean rotational velocity and squared velocity dis- 
persion are: 

(Vj) 2 = 2[< v\ >]edo ; (13) 

a 2 = a\ x + a\ 2 + cr 2 3 = 2[< a\ >] cdo + (1 - S)[< a\ >] cdo 

= (3-(J)[<«T?>]edo ; (14) 
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in terms of edge-on mean rotational velocity and velocity dispersion, and: 

2 

(v^f = — — [< v\ >] obs ; (15) 
sin i 11 

° 2 = 1 3 7 2 . [<gg>U ; ( 16 ) 

l — o COS z 11 

in terms of observed mean rotational velocity and velocity dispersion. 

The intrinsic rotation parameter, by use of Eqs. ( THl) . (fl5T) . ( |T6l) and f |45l) . 

/ 2^ _ (W) 2 _ 2 l-0~COsM<£g>]obs 

UJint " cr 2 "sin 2 * 3-5 [<<x 2 >] obs ' [U) 

and sample objects may be represented on the (Oexl) plane using Eqs. (TlOT) 
and (ED. 

The dimensionless energy: 

~~ GM ' (18) 

may be considered as the extent to which gravitation is balanced by cen- 
trifugal forces at the top major axis, ai, as —Fc/Fq ~ k. With regard 
to [Mio kpc Gyr] units, Mio = lO lo m , the constant of gravitation is G = 
4.493 10 4 M^o 1 kpc 3 Gyr" 2 and assuming as typical values M = 10 2 Mio, &i — 
lOkpc, a\ x = 4.493 10 4 kpc Gyr" 1 , Eq. (ITHj) yields k = 0.1 as expected for 
elliptical galaxies. On the other hand, observed quantities are related to the 
effective radius, R e , instead of the major semiaxis, a±, and the dimensionless 
energy has to be approximated as: 

[< a\ >} cdo R e 

« edo = m ; (19) 



1 e 



in connection with edge-on configurations. 

For the sample of elliptical galaxies under consideration (N = 16), the 
values of the following parameters are listed in Tab.[TJ the effective (half- 
light) radius, R e (CV08); the intrinsic mean equatorial tangential velocity 
component, Eq. ffT5"]) ; the intrinsic velocity dispersion, a, Eq. ffTB"]) ; the 
galaxy stellar mass within the effective radius, M e = M(R e ) = M/2, under 
the assumption that luminosity traces the mass (CV08); the dimensionless 
energy, K edo , Eq. (fl9l) : the inclination angle, i, of the best fitting two-integral 
Jeans model (SIV); the anisotropy parameter, 5, determined from the solu- 
tion of the dynamic models, supposed to be axisymmetric (SX); the intrinsic 
ellipticity, e, deduced from the computed inclination, Eq. (fTOi) . under the as- 
sumption of axisymmetric configurations (SX); the sample object rotation 
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Table 1: Parameters calculated from the data related to a sample (N = 16) 
of elliptical galaxies, extracted from larger samples of early-type galaxies 
investigated within the SAURON project (SIV, N = 25; SX, iV = 48). 
Values of dimensional quantities are expressed in [M 10 kpc Gyr] units (1 
kpc/Gyr=0.978 46 km/s; 1 km/s=l. 02201 kpc/Gyr), and angles in degrees. 
Column captions: (1) NGC number; (2) effective (half-light) radius, R e 
(CV08); (3) intrinsic mean equatorial tangential velocity component, 
Eq. (fT5i) ; (4) intrinsic velocity dispersion, cr, Eq. (|T6l) ; (5) galaxy stellar mass 
within the effective radius, M e = M(R e ) = M/2, under the assumption that 
luminosity traces the mass (CV08); (6) dimensionless energy, K edo , Eq. ffT9l ; 
(7) inclination angle, i, of the best fitting two- integral Jeans model (S IV); (8) 
anisotropy parameter, S, determined from the solution of the dynamic mod- 
els, supposed to be axisymmetric (SX); (9) intrinsic ellipticity, e, deduced 
from the computed inclination, Eq. (|T0|) . under the assumption of axisym- 
metric configurations (SX); (10) sample object rotation parameter, (xl)int, 
Eq. (1TTI) ; (11) adjoint configuration rotation parameter, xl, Eq. ©; (12) kine- 
matic classification, where F and S denote fast and slow rotators, respectively 
(S X). For the original data refer to the parent papers (S IV; S X). For further 
details refer, in addition, to earlier attempts (B05; CV08; C09). 
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parameter, (xDint, Eq. ([T7]) : the adjoint configuration rotation parameter, 
Xvi Eq- ©; together with the kinematic classification, where F and S denote 
fast and slow rotators, respectively (SX). For the original data refer to the 
parent papers (SIV; SX). For further details refer, in addition, to earlier 
attempts (B05; CV08; C09). It can also be seen that K edo = 0.20-0.34 for the 
whole sample, which implies gravitation is not balanced by centrifugal force 
in elliptical galaxies, as expected. 

The position of fast and slow rotators on the (Oexl) plane, is shown in 
Fig. [2] as squares and diamonds, respectively. Three curves are also repro- 
duced from Fig.fU namely equilibrium figures (£ = 1, full line) and nonequi- 
librium figures (£ = 1.1, up on the first quadrant; ( = 0.5, down on the first 
quadrant; dotted lines) with isotropic stress tensor. The caption of the ver- 
tical lines is the same as in Fig.[TJ The adjoint configurations with isotropic 
stress tensor and equal shape, are also positioned on the sequence of equi- 
librium figures using the same symbols, as the result of a vertical shift. A 
similar procedure holds, in general, for an assigned sequence of nonequilib- 
rium figures with isotropic stress tensor. An inspection of Fig. [2] shows the 
following features. 

All sample objects lie below the lower sequence (£ = 0.5), with the ex- 
ception of the most flattened galaxy (NGC 2974), which is placed just above 
the middle sequence (£ = 1). 

In general, slow rotators (diamonds) are placed on the left with respect 
to fast rotators (squares) with a single exception (NGC 3379) which, on the 
other hand, shows a nearly isotropic stress tensor (8 = 0.03), a nearly spher- 
ical shape (e = 0.04), and a velocity ratio, {[< vj >] bs/[< crjj >]obs}^ 2 = 
0.14, equal to the maximum value found in slow rotators (NGC 5813). A 
similar trend is exhibited by adjoint configurations with isotropic stress ten- 
sor, which are vertically shifted on the sequence of equilibrium figures, full 
line in Fig.[2j 

If slow rotators (including NGC 3379) are assumed to be nonrotating 
(v$ = 0) within the observational errors and elongated (e < 0) due to an 
anisotropic stress tensor (an = a 22 < 033 or 5 < 0), the related adjoint 
configurations with isotropic stress tensor are elongated due to imaginary 
rotation, conformly to Eq. More specifically, the major axis coincides 
with the polar axis and the meridional axis ratio exceeds unity, yielding a 
negative ellipticity, e = 1 — £31 < 0, which produces, in turn, a negative 
anisotropy parameter, 5 = 1 — crf 3 /cr n < 0, by generalizing Eq. (!7aj) to the 
case under discussion, o\\ = a 2 2 < ^33 • 

In this view, let efl a , 5^, be the intrinsic ellipticity from Eq. (1TU1) and the 
anisotropy parameter determined by comparison with the dynamic models, 
for axisymmetric flattened configurations (SX), and e e io, 8 e \o, their counter- 
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Table 2: The intrinsic ellipticity, e, the anisotropy parameter, 5, and the 
square root of the absolute value of the rotation parameter, Ix^p 2 , f° r sl° w 
rotators (NGC 3379 included) supposed to be flattened (fla) as in Tab.[U 
or elongated (elo) with equal major axis, according to Eqs. (120]) and (I2"T1) . 
The square root of the absolute value of the rotation parameter has been 
listed for better appreciating the difference between flattened and elongated 
configurations. 
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parts related to equal values of intrinsic axes and stress tensor components, 
but with the major axis related to the symmetry axis instead of the equato- 
rial plane. Accordingly, the following changes must be performed: £31 — > ; 
°33 "~ > Oil i Oil — * 033; which, using Eqs. (J7J) and ffTUj) . yields: 

e e io = 1 - eii 1 = 1 " (1 " efla)" 1 = ; (20) 

4.o = 1 - (^f 2 = 1 - (1 - M- 1 = -^t ; (21) 

Wll/ l-dfl a 

where the values related to flattened configurations are listed in Tab.[T] and 
again in Tab. [2] for slow rotators (NGC 3379 included), together with their 
counterparts related to elongated configurations and, in both cases, the values 
of the rotation parameter, Xv, for adjoint configuration with isotropic stress 
tensor. The square root of the rotation parameter has been listed for better 
appreciating the difference between flattened and elongated configurations. 

The position of slow rotators (NGC 3379 included) on the (Oexl) plane 
in the case under discussion, is indicated in Fig.[H] as diamonds and a single 
square, while nothing changes for the remaining fast rotators (squares) and 
captions, with respect to Fig. [21 
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4 Comparison with earlier attempts 



For equilibrium figures (£ = 1) viewed edge-on (i = 90°), the combination of 
Eqs. flS}, P, (El), and (ED yields: 

which is formally coincident with classical results, where <v^> and < > 
are intended as a mass-averaged, one-dimension squared rotation velocity 
and related variance, respectively (e.g., B05). 

If, on the other hand, < ujj > and < cr| > are conceived as sky-averages, 
according to the current notation, and the kinetic energy is decomposed 
into contributions from ordered and random motions, the right-hand side of 
Eq. (1221) has to be divided by a correction factor. The result is (B05): 

, [< Vj >]edo (1-^(^)11/(^)33-1 ,„,x 



where ip is an integral which depends on the density profile and the rotation 
curve of the matter distribution. 

In the present investigation, Eq. (1221) holds provided the kinetic energy 
is decomposed into contributions from cylindrical rotation and residual mo- 
tions, which implies a different anisotropy parameter with respect to Eq. (123H . 
The left-hand sides of Eqs. (12 2 j) and 



(Xw)edo — r 20 ' (^) 

[< 0| >Jcdo 

must necessarily coincide, if deduced from observations, which is true also for 
the shape factors. Accordingly, the combination of Eqs. (122 1) and ( J23l) yields: 



, (Bscl)ll (1 - ^b)(-Bsci)ii/(-Bsci)33 ~ 1 /n.x 

1 Cj (i?sel) 3 3 ^(l-5B)(i?scl)ll/(i?sel)33 + l ' ^ ' 

where the indices, C and B, denote anisotropy parameters related to the 
current and earlier (B05) attempt, respectively. 

The values of the anisotropy parameter for sample objects, listed in Tab.[H 
are taken from an earlier attempt (SX), where they have been determined in 
accordance with Eq. ( 1231) . For this reason, a comparison is needed between 
values of the rotation parameter, (xl)c and (xI)b, defined by Eqs. (l22i) and 
), respectively. To this aim, the rotation parameter ratio: 

xIb = ^ = HI-S b ) { 7 ^ + 1 ; 0<6 31 <1; (26) 



seU 33 
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shall be considered, which is maximum for isotropic stress tensor, 5 = 0. In 
the spherical limit, e 3i = 1, (-B se i)n = (-B sel ) 33 , which implies 5 = 0, and 
Eq. ( l26l) reduces to: 

Xcb = ^ + i ; e 3 i = i ; (27) 

in the flat limit, e 31 = 0, (B scl ) n = tt/2, (-B sel ) 33 = 0, which implies Xcb 
+oo unless 5 = 1. In dealing with elliptical galaxies, e 3 i > 0.3, the flat limit 
is never attained. For realistic density profiles, numerical computations yield 
ip m 0.131 in presence of flat rotation curves, v$ = const (B05), and ip m 0.15 
using Jeans models (SX). 

The results are listed in Tab. El where the values of the following pa- 
rameters are calculated or deduced from Tab.Q] for sample objects: the 
intrinsic meridional axis ratio, e 3 i = 1 — e, Eq. (JTUl) : the shape factor ra- 
tio, 623 = (-Bsei)22/(-B se i)33, Eq. (l3bl : the inclination angle, i, of the best 
fitting two- integral Jeans model (SIV); the anisotropy parameter, 5, deter- 
mined from the solution of the dynamic models, supposed to be axisymmetric 
(SX); the sample object rotation parameter, (xl)cdo = [< vf\ >]cdo/[< cr?\ > 
]edo, Eqs. (JTTJ and (TT2T) : the adjoint configuration rotation parameter, (xl)c, 
Eq. (1221) : the adjoint configuration rotation parameter, (xI)b, Eq. (12B1 ; the 
rotation parameter ratio, Xcb> Eq. f }2"6"]) ; the ratio, p = 1/(\/2xcb); together 
with the kinematic classification, where F and S denote fast and slow rota- 
tors, respectively (SX). For sample objects with meridional axis ratio below 
the bifurcation point, e 3 i < (e 3 i)bif = 0.582 724, the shape factors are deter- 
mined for triaxial instead of axisymmetric configurations, where the minor 
equatorial axis is assumed to be parallel to the line of sight. 

An inspection of Tab. [H] shows that Xcb =1-15-1.23 and 1/(v^2xcb) = 
0.64-0.66. Curiously, the last range is consistent with the squared coefficient, 
2 x 0.57 2 , of the best fitting relation between rotation parameters derived from 
integral- field and long-slit stellar kinematics (SX), the second expressed in 
terms of one-dimension peak velocity. 

5 Discussion 

Though the generalization of the ellipticity-rotation plane, (Oe^), takes 
nonequilibrium figures into consideration, still galaxies may be thought of as 
fully virialized (( = 1) unless a major merger is going on. On the other hand, 
cluster of galaxies are presently assembling and radial motions dominate at 
large clustercentric distances, which implies dynamical evolution (C !)■ F° r 
this reason, only the sequence of equilibrium figures has been considered in 
dealing with elliptical galaxies. The rotation parameter, Xvi depends on the 
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Table 3: Comparison between intrinsic rotation parameters, (xl)c and (xI)b, 
calculated for sample objects listed in Tab.HJ using Eqs. (1221) and (123"]) . re- 
spectively, for isotropic stress tensor (5 = 0) and if) = 0.15 (SX). Also 
listed are values of the ratios, Xcb = (xl)c/ (xI)b and p = 1/(v2~Xcb)- 
The remaining parameters are taken or deduced from Tab.HJ For merid- 
ional axis ratios below the bifurcation point, e 3 i < (e3i)bif = 0.582 724, the 
shape factors are determined for triaxial instead of axisymmetric configura- 
tions, where the minor equatorial axis is assumed to be parallel to the line of 
sight. Column captions: (1) NGC number; (2) intrinsic meridional axis ratio, 
e 3 i = 1 - e, Eq. (JTHJ); (3) shape factor ratio, 623 = (^1)22/(^1)33, Eq. ((30); 
(4) inclination angle, i, of the best fitting two- integral Jeans model (S IV); (5) 
anisotropy parameter, S, determined from the solution of the dynamic mod- 
els, supposed to be axisymmetric (S X); (6) sample object rotation parameter, 
(xl)edo = [< v\ >]cdo/[< of >]edo, Eqs. §H§ and CE2J); (7) adjoint configura- 
tion rotation parameter, (xl)c, Eq. ( 1221) ; (8) adjoint configuration rotation 
parameter, (xI)b, Eq. (123]) ; (9) rotation parameter ratio, Xcb> Eq. (126]) ; (10) 
ratio, p = 1/(v2xcb); (H) kinematic classification, where F and S denote 
fast and slow rotators, respectively (SX). 
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intrinsic squared mean equatorial tangential velocity component and velocity 
component dispersions, according to Eq. (jSJ). 

A different rotation parameter has been defined in an earlier attempt 
(B05), which depends on the sky-averaged mean streaming velocity parallel 
to the line of sight, the sky- averaged streaming velocity (parallel to the line 
of sight) dispersion, and the sky-averaged component velocity (parallel to the 
line of sight) dispersion related to nonstreaming motions. 

The inclination angle, i, and the anisotropy parameter, 5, have been de- 
termined for sample objects by comparison with the dynamic models, under 
the assumption of axisymmetric configurations (oi = a^), axisymmetric stress 
tensor (an = (T22), and flattened shapes (ai > a 3 , a n > a 33 ), as outlined in 
the parent paper (SX). On the other hand, it is shown in Fig. [5] that four 
sample objects are predicted to be triaxial as the ellipticity value related to 
the bifurcation point (ebif = 0.417276), marked by the vertical dotted line, 
is exceeded, and the most flattened configurations could be barlike (e.g., 
e 2 i = 0.432 232 for e 3 i = 0.345 069). Then some caution must be adopted in 
the interpretation of the trixiality, expecially for the most aspherical sample 
object (NGC 2974). To this respect, it is worth remembering that galax- 
ies are embedded in dark (nonbaryonic) matter haloes, according to current 
cosmological scenarios. 

The presence of a massive, embedding subsystem stabilizes the inner 
spheroid and shifts the bifurcation point (from axisymmetric to triaxial 
configurations) towards increasing ellipticities (Durisen 1978; Pacheco et al. 
1986). In the special case of homogeneous subsystems, the bifurcation point 
is attained at e ~ 0.7 for comparable masses within the volume of the inner 
spheroid (Caimmi 1996a), which can be considered as an upper limit. In 
fact, by comparison with the dynamic models, fast rotators appear to be ax- 
isymmetric while, paradoxically, slow rotators exhibit a moderate triaxiality 
(SX). Using triaxial dynamic models could provide better understanding on 
this point. An inspection of Fig.[2]and Tab. CD shows that, for slow rotators 
(including NGC 3379) < e < 0.2; < 5 < 0.15; < xl < 0.15; and 
for fast rotators (excluding NGC 3379) 0.2 < e < 0.65; 0.15 < 5 < 0.35; 
0.15 < xl < 0-5; which could be an alternative kinematic classification with 
respect to earlier attempts (e.g., Emsellem et al. 2007). Richer samples 
should be dealt with to provide more conclusive evidence on this point. 

If slow rotators (including NGC 3379) are nonrotating within the obser- 
vational errors, and their shapes are elongated due to negative anisotropy 
parameters, 5 < i.e. an = a 2 2 < 033, a different kinematic classification 
is shown in Fig. [3] and Tab. [21 With regard to adjoint configurations with 
isotropic stress tensor, fast and slow rotators are affected by real and imag- 
inary systematic rotation, respectively. It is, of course, a limiting situation, 
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in the sense that no zone of avoidance is expected on the sequence of adjoint 
configurations. If nonrotating elongated configurations really exist, slow ro- 
tators should be divided into two subclasses, namely (1) flattened in real 
rotatin, and (2) elongated in imaginary rotation. In the latter alternative, 
elongated instead of flattened dynamic models should be used for determin- 
ing the inclination angle and the anisotropy parameter by comparison with 
the data (S IV; S X), to gain consistency and to test the above interpretation. 

Given a spherical galaxy with isotropic stress tensor, rotation (around 
a selected axis) kinetic energy may be added in infinite ways between two 
limiting situations, namely (a) systematic rotation, where circular motions 
are either clockwise or counterclockwise, and (b) random rotation, where the 
mean circular velocity is null. In any case, all the observables remain un- 
changed with the exception of the angular momentum and the stress tensor. 
More specifically, any system may be related to its adjoint configuration with 
isotropic stress tensor, without loss of generality, for ellipticity values above 
a threshold (Appendix [B]) . 

In this view, the elliptical side of the Hubble sequence may be inter- 
preted as a sequence of equilibrium (adjoint) figures where the ellipticity, e, 
is increasing with the rotation parameter, xl- The original (elliptical side of 
the) Hubble sequence may be generalized on two respects, namely (1) from 
equilibrium (£ = 1) to nonequilibrium (£ ^ 1) figures, and (2) from real 
(e > 0, — 0) to imaginary (e < 0, xl < 0) rotation. More specifically, 
real rotation is related to flattened configurations and imaginary rotation 
to elongated configurations, spinning around the minor and the major axis, 
respectively. 

The above classification is grounded on the assumption of homeoidally 
striated density profiles, regardless of the mechanism of formation. More 
specifically, configurations of equal shape related to different assembling pro- 
cesses belong to the same class. In reality, isophotes in elliptical galaxies may 
be boxy i.e. "overelliptic" or disky i.e. "underelliptic" where, in any case, 
the axis ratio changes with radius. This dichotomy, together with the one 
related to dwarf ang giant elliptical galaxies, is interpreted as due to different 
past histories (e.g., Kormendy et al. 2009). On the other hand, a simplified 
description in terms of isophotes with constant axis ratio at all radii, implies 
a classification of equilibrium (or nonequilibrium with fixed virial index) fig- 
ures which is independent of the formation and evolution processes. 
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6 Conclusion 



The results of earlier investigations on homeoidally striated MacLaurin spheroids 
and Jacobi ellipsoids (CM05; Caimmi 2006a; C07) are used in the current 
attempt for the representation of nonequilibrium figures on the ellipticity- 
velocity plane, and of figures in imaginary rotation, where the effect is elon- 
gating instead of flattening, with respect to the rotation axis. The key con- 
cept is that the addition of kinetic energy related to tangential equatorial 
velocity components makes distorted boundaries regardless of what fraction 
translates in an increment of angular momentum, and what fraction in an in- 
crement of stress tensor equatorial components. Then any system admits an 
adjoint configuration with isotropic stress tensor, for ellipticity values above a 
threshold (Appendix [B]), and the related sequence can be considered without 
loss of generality. The kinetic energy is decomposed into contributions from 
cylindrical rotation, E vot = M(v^) 2 /2, and tangential equatorial component 
velocity dispersion plus residual motions, E res = M^a^ + o^ w + <r| 3 )/2, and 
the rotation parameter is defined as xl — E I0t /E Ie8 . 

Elliptical galaxies are idealized as homeoidally striated MacLaurin spheroids 
and Jacobi ellipsoids, and their position on the ellipticity- velocity plane is 
inferred from observations related to a sample (CV08, N = 16) extracted 
from richer samples of early-type galaxies investigated within the SAURON 
project (SIV, = 25; SX, N = 48). The location of model galaxies on 
the (Oexl) plane is determined through the following steps (C09): (i) se- 
lect SAURON data of interest; (ii) calculate the parameters appearing in 
the virial equations; (iii) make a correspondence between model galaxies and 
sample objects; (iv) represent model galaxies on the (Oexl) plane. 

The main results found in the present investigation may be summarized 
as follows. 

(1) Sequences of homeoidally striated MacLaurin spheroids and Jacobi ellip- 

soids with isotropic stress tensor are defined and plotted on the (Oexl) 
plane, without loss of generality: for any system with anisotropic stress 
tensor an adjoint configuration exists, with isotropic stress tensor and 
remaining parameters unchanged with the exception of the angular 
momentum, for ellipticity values above a threshold. 

(2) Sequences of homeoidally striated MacLaurin spheroids and Jacobi el- 

lipsoids are generalized to nonequilibrium figures. 

(3) Sequences of homeoidally striated MacLaurin spheroids and Jacobi el- 

lipsoids are generalized to imaginary rotation which acts in elongating, 
instead of flattening, with respect to the rotation axis. 
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(4) An alternative kinematic classification with respect to earlier attempts 

(e.g., Emsellem et al. 2007) is proposed, where slow rotators are char- 
acterized by low ellipticities (0 < e < 0.2), low anisotropy parameters 
(0 < 5 < 0.15), and low rotation parameters (0 < xl < 0.15), and 
fast rotators by large ellipticities (0.2 < e < 0.65), large anisotropy 
parameters (0.15 < 5 < 0.35), and large rotation parameters (0.15 < 
xl < 0.5). Richer samples should be used to test the validity of the 
above interpretation. 

(5) A possible interpretation of slow rotators as nonrotating and elongated 

due to a negative anisotropy parameter, instead of flattened due to a 
positive anisotropy parameter, is exploited. 

(6) The elliptical side of the Hubble morphological sequence is interpreted 

as a sequence of equilibrium (adjoint) figures where the ellipticity is an 
increasing function of the rotation parameter. Accordingly, slow rota- 
tors correspond to early classes (E0-E2 in the oblate limit and E— 2-EO 
in the prolate limit) and fast rotators to late classes (E3-E6). In this 
view, boundaries are rotationally distorted regardless of what fraction 
of tangential equatorial velocity components is related to angular mo- 
mentum, and what fraction to stress tensor equatorial components. 
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Appendix 

A Rotation and residual kinetic energy 

With regard to a generic matter distribution, let (0x1X2X3) be a reference 
frame where the origin coincides with the centre of mass, and the coordi- 
nate axis, X3, coincides with the rotation axis. Let v<f,(xi, X2, X3, t) be the 
tangential equatorial (with respect to the rotation axis) velocity component, 
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related to the generic particle at the point, P(x±, x 2 , x 3 ), at the time, t. The 
corresponding angular velocity is defined by the relation: 

v (f ,(x 1 ,x 2 ,x 3 ,t) = Q(x 1 ,x 2 ,x 3 ,t)w ; w 2 = x\ + x\ ; (28) 

which can also be extended to mean values, v^(xi,X2,x 3 ,t), Q(xi,x 2 ,x 3 ,t), 
whithin an infinitesimal volume element, dxi dx2 dx 3 , placed at the same 
point, P. In the following relations, the time dependence shall be omitted on 
the quantities at the left-hand side, to gain simplicity. 

The mass-weighted tangential equatorial velocity component and squared 
tangential equatorial velocity component read: 

^ = ij/f// /^( Xl ' X2 ' :r 3,i)p(a;i,^2,a;3,i)da;ida;2da;3 ; (29) 

^ = ~M 1 1 J v l( x ^^ x ^^ x 3^ t )p( x i^ x 2^ x 3,t)dx 1 dx 2 dx 3 ; (30) 

where M is the total mass and p the local density. The empirical variance 
is: 



°U = (4) ~ (w) 2 ; (si) 

by definition. 

The moment-of-inertia-weighted angular velocity and squared angular ve- 
locity read: 

H= | / f /n(x 1 ,x 2 ,,3,^Vfe,^,X3, i )d, 1 dx 2 dx3 ; (32) 
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(VP) = — J J J n 2 (xi,x 2 ,x 3 ,t)w p(xi,x 2 ,x 3 ,t) dxi dx 2 dx 3 ; (33) 

where /3 is the moment of inertia with respect to the rotation axis, x 3 . The 
empirical variance is: 

° 2 nn =W)-(W ; (34) 

by definition. 

A link between mass and moment of inertia is provided by the following 
relation: 

*-#-i///^.W)d.^ i ,35) 

where Rq 3 is the radius of gyration with respect to the rotation axis, x 3 . 
Using Eqs. (J2HD and (JHSD, the combination of Eqs. (|3TJ|), (l3T]h (133]), and 
yields: 



(«?) = (^G3^ (36) 
(vj) 2 + ols = [(W + a 2 nn }R 2 G3 ; (37) 
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while, on the other hand, v$ ^ QRg3, a <t><t> a no,RG3i an d a different averag- 
ing is needed on either or Q. 

To this aim, the (Mi?Q3)-weighted circular velocity, vu, and the (MRq 3 )- 
weighted angular velocity, fl v , are defined as: 

vn = QRG3 = ^r-f I [ v^(xi,x 2 ,x 3 ,t)wp(x 1 ,x 2 ,x 3 ,t)dx 1 dx 2 dx 3 ; (38) 



I. 



3 



Vl v = = / [ [ {l(x 1 ,x 2 ,x 3 ,t)wp(x 1 ,x 2 ,x 3 ,t)dx 1 dx 2 dx 3 ; (39) 

-r<G3 1 3 J J J 

and the angular momentum with respect to the rotation axis, x 3 , is defined 
by the integral on the right-hand side of Eq. ( J32l) or (l38l) , which yields: 

J 3 = m = MvUR G3 ; (40) 

in terms of the moment-of-inertia- weighted angular velocity and the (MRq 3 )- 
weighted circular velocity, respectively. 

The combination of Eqs. ( 1351) . ( 1361) . ( 1381) . and ( 1391) yields the following 
expression for the radius of gyration: 



R 2 h = ^ == {^ = ^_ . (41) 



where the differences, (v^) — (vq) 2 , (Q 2 ) — (Q v ) 2 , cannot be related to empirical 
variances as the corresponding terms are weighted in different ways. 

Accordingly, the rotation kinetic energy related to tangential equatorial 
velocity components may be expressed as: 
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(£ k inW = 2 M [W + a U\ = + 2 M °U ; ( 42a ) 

E rot = ^M(^) 2 ; (42b) 

and the total kinetic energy reads: 

-^kin = E rot + E ies ; (43a) 

; (43b) 



£res = ~M 



K) + (vi) + al 



{yj) = M + °ww = °ww ; ( 43c ) 

[vl) = {vif + al 3 = al 3 ; (43d) 

where v w and v 3 are the radial equatorial and polar velocity components and 
= 0$ = provided the centre of mass coincides with the origin of the 
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coordinates (for further details refer to e.g., C07). Accordingly, Eq. f!43bj) 
reduces to: 

E ves = l -Ma 2 = l -M{a^ + al w + al i ) ; (44) 
which is equivalent to: 

^ = ^ = -¥(4 + 4 + 4) ; (45) 

in Cartesian coordinates. 

B Systematic and random rotation excess 

Let two matter distributions be characterized by equal density profiles and 
shapes, but isotropic and anisotropic stress tensor, respectively, and different 
amount of angular momentum. The generalized virial equations, Eqs. (JT]), 
read: 

2[(-E r ot) PP ]ani ~\~ 2 [(-£/res)pp]ani 

+ (£ P ot) PP = ; p= 1,2,3 ; (46a) 
2[(£ rot ) PP ]iso + 2[(£ rcs ) pp ] iso + (E pot ) pp = ; p = 1,2,3 ; (46b) 
(E Iot ) n + (E vot ) 22 = E Iot ; (£ rot ) 33 = ; (46c) 

where (E res ) pp = ( pp E res is the residual kinetic-energy tensor of the related 
equilibrium figure. 

In both cases, the potential-energy tensor remains unchanged, which im- 
plies the following relation: 

[(-E'res)33]ani = [(-E'res) PP ]iso ', P = l,2,3 ] (47) 

and the substitution of Eq. (F4T1) into (14*61 yields: 

qq\\so 

[(-^res)gg]iso + A(E res ) qq 
= [(E vo t) qq ] iso +[(E vcs ) qq ] iso ; g = l,2 ; (48a) 

A(E Tot ) qq = [(-E'rot)q 9 ]ani ~ [(E T0 t) qq\iso ] (48b) 
^(-E'res)gg [(-^res)gg]ani 

- [{E res ) qq \ iso ; (48c) 

which is equivalent to: 

A(E TOt ) qq + A(E Tes ) qq = ; (49) 

where the first and the second term of the sum may be conceived as a (positive 
or negative) systematic and random rotation excess, respectively. Random 
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rotation has to be intended as related to a selected axis, with a null mean 
value. 

With regard to nonequilibrium figures, the validity of the relation: 

A(E TOt ) qq + A(E rea ) qq = ; (50) 

means that a fixed amount of rotation energy has been converted into residual 
energy (or vice versa). The combination of Eqs. ( 1491) and (1501 yields: 

A(E ies ) qq = A(E ies ) qq ; (51) 

and Eq. (1501) may be cast under the explicit form: 

A[(v^) qg } 2 + A(a 2 J=0 ; (52) 

according to Eqs. fpf2l and (HHI) . The above results may be reduced to a single 
statement. 

Theorem Given a matter distribution with assigned density profile, shape, 
virial index, ( = E res /E res , and isotropic stress tensor, an infinity of ad- 
joint configurations with anisotropic stress tensor exist, for which the sum of 
systematic and random rotation excess is null. 

According to earlier attempts (Lynden-Bell 1960, 1962; Meza 2002) the 
distribution function is independent of the sign of the tangential velocity com- 
ponents, and the whole set of possible configurations is characterized by an 
equal amount of both kinetic and potential energy. In other words, clockwise 
and counterclockwise circular motions are indistinguishable to this respect, 
passing from systematic (maximal squared mean tangential equatorial ve- 
locity component) to random (minimal squared mean tangential equatorial 
velocity component) rotation. An anisotropic stress tensor could be due to 
the prevalence of either systematic rotation (an = 022 < 033) or random 
rotation (an = 022 > 033) with respect to the adjoint configuration with 
isotropic stress tensor. 

In the limiting case of flat configurations, e 3i — > 0, (-E po t)33 — * (e.g., 
Caimmi 2009a), Eq. (Tib]) implies (E TCS ) 33 = C33-^res -> and (E res ) qq -> 
for isotropic stress tensor. In other words, flat configurations with isotropic 
stress tensor must necessarily be self-gravitating. Rotation curves depend on 
density profiles, ranging from linear curves related to homogeneous systems 
to Keplerian curves related to Roche systems. Accordingly, a^ = in the 
former alternative, growing up to a maximum, a^ > 0, in the latter. But 

> would imply (E res ) qq > 0, which is in contradiction with either suffi- 
ciently flattened configurations or isotropic stress tensors. Then a threshold 
in meridional axis ratio exists, €31 = e^, below which configurations with 
isotropic stress tensor cannot exist for an assigned density profile. 
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Figure 1: The rotation paramerer, Xvj as a function of the meridional el- 
lipticity, e = 1 — £31, for sequences of equilibrium (full) and nonequilib- 
rium (dotted) figures with isotropic stress tensor. For equilibrium figures, 
C = E res /E res = 1. For nonequilibrium figures, ( = 1.1, 1.2, 2.0 (above the 
equilibrium sequence on the first quadrant), and ( = 1/1.1, 1/1.2, 1/2.0 
(below the equilibrium sequence on the first quadrant). Configurations ly- 
ing on the first and third quadrant are in real and imaginary rotation, re- 
spectively. Nonrotating configurations are placed at the origin. Caption of 
vertical lines: dotted - bifurcation from axisymmetric to triaxial configura- 
tions; dashed - bifurcation from triaxial (left) and axisymmetric (right) to 
pear-shaped configurations; dot-dashed - onset of dynamical instability in 
axisymmetric configurations. 
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Figure 2: Position of elliptical galaxies listed in Tab.HJ on the (Oexl) plane, 
denoted as squares and diamonds for fast and slow rotators, respectively. 
Curves are taken from Fig. CD, related to equilibrium figures (£ = 1, full line) 
and nonequilibrium figures (£ = 1.1, up on the first quadrant; ( = 0.5, down 
on the first quadrant; dotted lines) with isotropic stress tensor. Caption of 
vertical lines as in Fig.[TJ The adjoint configurations with isotropic stress 
tensor and equal shape, are also positioned on the sequence of equilibrium 
figures using the same symbols, as the result of a vertical shift. A similar pro- 
cedure holds, in general, for an assigned sequence of nonequilibrium figures 
with isotropic stress tensor. 
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Figure 3: Same as in Fig. El but under the assumption that slow rotators 
(NGC 3379 included) are nonrotating within the observational errors and 
elongated due to an anisotropic stress tensor. For further details refer to the 
text. Captions as in Fig. El 
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